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A theory for the microscopic structure and the vibrational properties of soft sphere glass at finite temperature 
is presented. With an effective potential, derived here, the phase diagram and vibrational properties are 
worked out around the Maxwell critical point at zero temperature T and pressure p. Variational arguments and 
effective medium theory identically predict a non-trivial temperature scale T* ~ p( 2 ~ a )/U- a ) with a « 0.17 
such that low-energy vibrational properties are hard-sphere like for T > T*, and zero-temperature soft-sphere 
like otherwise. However, due to crossovers in the equation of state relating T, p , and the packing fraction </>, 
these two regimes lead to four regions where scaling behaviors differ when expressed in terms of T and </>. 
Scaling predictions are presented for the mean-squared displacement, characteristic frequency, shear modulus, 
and characteristic elastic length in all regions of the phase diagram. 


A wide range of amorphous materials including granu¬ 
lar materials, foams, molecular glasses and colloids ex¬ 
hibit a transition from liquid-like to solid-like behav¬ 
ior. In the solid phase, these materials display anoma¬ 
lous elastic properties. In particular, amorphous solids 
universally present an excess of vibrational modes over 
the Debye model (that treats vibrational modes as plane 
waves) a phenomenon referred to as the ‘boson peak 11 . 
The boson peak affects various properties of the solid, 
including heat transport 2 Si as well as the spatial hetero¬ 
geneity of the elastic response 7 AS). Its presence indicates 
that glasses lie close to an elastic instability, suggest¬ 
ing a connection with the glass transition where rigidity 
emerged 11 l 12 l For these reasons this phenomenon is in¬ 
tensely studied^. _ 

It was proposed that in a variety of glassed ^ the bo¬ 
son peak is controlled by some aspects of the short-range 
microscopic structure, in particular a measure of the 
particle connect ednes4^““ as well as t he cha racteristic 
force with which nearby particles interaei P 9 l 2Q l both well- 
known to affect the stability of engineering structured^ 1 . 
Detailed predictions on transport and sound dispersion 
can be made in this framework 15 , which have been re¬ 
cently supported by experiments 22 . In this approach, 
particles interacting via a purely repulsive potential at 
zero temperature^HSS are predicted to display singular 
vibrational properties: the spectrum has a characteris¬ 
tic frequency w* satisfying oj* ~ 2 — z c near the jam¬ 
ming transition 1 ^ as confirmed numerically®. Here 2 is 
coordination, measured by the average number of con¬ 
tacts per particle, and z c = 2d is the minimal coordi¬ 
nation for mechanical stability in dimension d, the so- 
called Maxwell threshold. The boson peak frequency 
u>bp, however, is essentially zero, with to bp "C w* in¬ 
dicating that the system is marginally stable 2 ^ and im¬ 
plying that z — z c ~ yj(j> — <j) c , where </> e is the jamming 
packing fraction where pressure vanishes (that can de¬ 
pend on the configuration considered). 


Since in practice interaction potentials are always an- 
harmonic, one m ay qu estion if these results survive at 
finite temperature 50 37 Hard spheres, arguably the sim¬ 
plest glass former, are a particularly challenging case for 
which the potential is discontinuous, implying that the 
Hessian of the energy (whose eigenvectors correspond to 
the normal modes) is not well-defined. However, it was 
argued that a coarse-grained free energy can be com¬ 
puted in that can 32 All jf 0 ne can find times r such that 
r c -C t <C r a , where r c is a microscopic collision time, 
and r a is the relaxation time of the structure (this con¬ 
dition is always achieved in the glass phase). This free 
energy captures the volume of phase space around a given 
meta-stable state (or “vibrational entropy”), and can al¬ 
ways be expressed in terms of the mean particle position 
within a meta-stable state. This approach leads to an 
effective interaction V e // between particles, which is in 
general multi-body. Near the jamming point, however, 
this interaction is simply two-body (as recalled below), 
and only occurs between particles interacting (i.e. collid¬ 
ing) on the intermediary time r, leading to the definition 
of a contact network and allowing to extend the notion of 
coordination to hard spheres. The Hessian of the coarse¬ 
grained free energy describes the linear response, in par¬ 
ticular how small external forces affect the mean particle 
positions. Its eigenvectors define normal modes whose 
thermal fluctuations are inversely proportional to their 
associated eigenvalue, as confirmed numerically 35 . 

Theoretical arguments using such an effective poten¬ 
tial were used to predict scaling properties of the spec¬ 
trum of the H essian , as well as spatial properties of 
the eigenvector d 33 * 34 l One finds that the normal modes 
slightly differ from the soft sphere case: in both cases 
modes are extended and heterogeneous, but for hard 
spheres soft modes tend to distort a small amount of 
bonds (which are weak), whereas for soft spheres modes 
tend to distort all the bonds in the system. We thus coin 
hard-sphere soft modes as sparse , whereas for soft spheres 
we keep the previously used terminology of anomalous 
modes. This distinction leads to a small but detectable®! 
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(T/M) 1 ' 2 

M K T\5(f)\- K 
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5(f) 

—5(f) 

5(f) 1 ' 2 

T\5(f>\- 1/2 | <50| 1/2 

5(f>- 1/4 

Soft-entropic (SE) 

T < T* 

5(f) > 0 

5(f) 

—5(f) 

5(f) 112 

(T/\5(f>\) K T 1 ~ K \5(f)\ K 

(T'/<50) (a—1 

Anharmonic (AH) 

T » 54> 2 


VT 

VT 

yi/4 

rp/c/2 

rp(a-l)/S 


TABLE I. Regimes of the soft-sphere glass. Shown is the scaling behavior of pressure p, typical time-averaged gap (h), 
characteristic frequency u j*, mean-squared displacement ( 5R 2 ), intra-state shear modulus p, and correlation length of the 
normal modes l s (u*) versus temperature T and volume fraction deviation 5(j) = <j> — (f > c . As discussed in the main text, the 
soft regime contains a transition at a nontrivial temperature scale T* ~ < 5 <(/ 5+3ee F( 2 + 2S =) ^ 5(f) 2 ' 2 , invisible in the effective 
potential and equation of state (and therefore p and (h)), but strongly affecting vibrational properties. The exponents k = 
(4 + 20 e )/(3 + 0 e ) « 1.41 and a = (1 — 9 e )/(3 + 6 e ) ~ 0.17, where 6 e « 0.42 characterizes the force distribution of “extended” 
contacts at jamming- 1 ®. 


difference in vibrational properties. For the coordina¬ 
tion, it was predicted that 2 — z c ~ (cj> c — (j>) 0A1 . Other 
predictions, such as for the mean-squared displacement, 
were found to ag ree w ith recent replica calculations in in¬ 
finite dimension j 27 F 8 l, supporting the validity of the two 
approaches. 

The asymmetry between soft and hard spheres appar¬ 
ent in the scaling of coordination is surprising, since one 
can go continuously from one case to the other by con¬ 
sidering soft spheres, and by allowing to vary both (f> and 
2)3 0131 ], 3 6HH] The limit </> > <f> c , T —> 0 corresponds to soft 
spheres whereas hard spheres are recovered in the limit 
(j) < </> c , T —► 0. Previous theoretical attempt^ 39142 * 43 1 to 
describe vibrational properties in the ((f), T) plane did not 
include this necessary asymmetry. Here we extend the 
effective potential approach to the case of thermal soft 
spheres, and use our real-space description of vibrational 
modes to describe the phase diagram entirely. Our main 
results are summarized in the phase diagram of FigJT] 
and Table [T| We find that low-energy vibrational prop¬ 
erties display a single crossover scale, T*, when written 
in terms of pressure p and temperature T ; this transition 
is indicated by the shaded color change from red to blue 
in FigJT] However, if pressure is eliminated in favor of 
volume fraction cf>, this introduces two more transitions, 
reflecting crossovers in the equation of state. In sum, this 
gives four scaling regimes in the phase diagram. 

The definition of a contact network and effective po¬ 
tential we provide complements the approach of Ref. ED, 
who discussed anharmonicity in amorphous solids from 
an experimentalist’s point of view. The experimental 
question is: given the displacement covariance matrix 
for an anharnroni c syst em, what can one learn about the 
vibrational mode^BD? The authors of Ref. [3T] defined 
the ‘shadow system,’ which is the Hamiltonian system of 
harmonic springs whose stiffnesses are related to displace¬ 
ments by the usual equilibrium formula. The obtained 
effective stiffnesses imply an interaction between particles 
on the averaging time-scale used to obtain displacement 
correlations. As we will see, we find a remarkable sim¬ 
ilarity between our effective potential and the inferred 
potential obtained from particle displacements in the nu¬ 
merical simulations of Ref. m 


I. EFFECTIVE POTENTIAL 

A. Derivation 

We consider a system of N soft spheres at constant 
pressure p that is fluctuating around a metastable state. 
Thermal averages taken around this state are written by 
(•), while averages over the contacts will be denoted by ~. 
We assume that the coordination of the network of inter¬ 
acting particles is z c = 2d, i.e the system is isostatic, as is 
the case at the jamming transition p, T -4 0. The compu¬ 
tation of the vibrational free energy in such a meta-stable 
state is detailed in Appendix A, and here we sketch the 
argument. Isostaticity implies that the number of de¬ 
grees of freedom in the particle displacements {Sri}, dN, 
is precisely equal to the number of contacts, Nc = zN/2. 
Thus the partition function of the nretastable state, orig¬ 
inally an integral over the displacements {5r)}, can be 
written as an integral over the gaps {h a } between con¬ 
tacting particles (with h a < 0 for overlap). Assuming 
that all |<5rj| <C 1, the map from {Sfi} to {h a } is linear. 

To compress or dilate the system from a volume V c 
(where all particles are just touching) to V requires work 
W = —p(V — V c ). Stability implies that the time- 



FIG. 1. Phase diagram for marginally-stable harmonic soft 
sphere systems near the Maxwell point at (T,(f>) = ( 0 ,</> c ). 
The dashed line indicates the scale T ~ |<5^f>| 2 , while the 
solid line indicates T* ~ j^,( 5 + 3 ®e)/(2+2e e )_ j n re gi mes Hard 
Spheres (HS), Anharmonic (AH), and Soft Entropic (SE), the 
softest modes are ‘sparse,’ as in hard sphere systems, but with 
different scalings for the pressure: p ~ T/\5(f>\ (HS), p ~ VT 
(AH), and p ~ 5(j) (SE). In regime SO, the softest modes are 
‘anomalous,’ like T = 0 soft solids. Note that the lines are 
schematic. 
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FIG. 2. Effective stiffness at indicated temperatures (a) 
as predicted from the simple form Eq. (3f, and (b) as com¬ 
puted from the displacement covariance matrix in numerical 
simulations 31 ^ at <f> — <f> c = 0.02. 


averaged forces f a between particles satisfy force balance; 
this can be written W = — fah a , where the sum is 
over all contacts a of the metastable state, and we have 
used the virtual work principle®. Since W and the elas¬ 
tic energy U are sums over contributions from different 
contacts (we consider only pa ir pot entials), this leads to 
a single-gap partition function 32 33 

Z(/3, f) = e~ pG = J dh e -mh) e -pfh, 

where the force / fixes the time-averaged gap ( h) 
through (h) = dG/df , and /3 = 1 /(hgT). The effec¬ 
tive potential 14// is obtained by a Legendre transform 
V e ff((h)) = G(f ) — f(h), from which it follows that / = 
—dV e ff((h))/d(h). We consider a finite-range harmonic 
potential U(h) = \e\h/a \ 2 when — \a < h < 0, and 0 
otherwise. We now take units in which e = cr = kg = 1. 
As described in Appendix A, it is simple to analytically 
extract the limiting behavior: when / <C VT, (h) « T //, 
whereas when / VT , ( h ) ~ — /. These correctly re¬ 
duce to the effective hard-sphere potential from Refs. (32] 
and[331 in the small-force limit, and the original harmonic 
potential in the large-force limit, respectively, and more¬ 
over they indicate the cross-over scales. In the interme¬ 
diate regime, |(/i)| < VT, there is a smooth cross-over 
between these behaviors. 

Since our goal is to obtain scaling behavior, it is conve¬ 
nient to have a simple form that incorporates the salient 
features of 14//, without adding extra structure. The sim¬ 
plest is the one satisfying ( h ) = —/ + T/f, which leads 
to an effective force law 

/(</i}) = 4</i} + ±y</i) 2 +4T, (2) 


introduction, it was argued that the effective potential 
can be measured by displacement covariances using the 
usual harmonic formula 31 . To test this claim explicitly, 
in Fig. [2b we reproduce the effective stiffness k e ff ex¬ 
tracted from displacement covariances in Ref. |3T from 
numerical simulations of a thermal soft-sphere system in 
d = 2. A very good qualitative agreement is obtained, 
with no fitting parameters. Moreover, in Ref. 31] it was 
observed that by varying 5(f) = <f> — (j) c from ss —0.02 to 
+0.05, the effective stiffness function k e ff((h)) remained 
constant, and the different systems merely sampled dif¬ 
ferent ranges of ( h ); this feature is reproduced by Eq. <[3j) , 
which has no explicit dependence on <f>. This observation 
supports that corrections to the effective potential away 
from isostaticity are small in this range of 5(f), as was 
checked previously for hard spheres 32 . 


B. Force Distribution 

The disordered geometry of the metastable state is 
characterized by the statistics of time-averaged particle 
positions, which determines in turn the distribution of 
forces between particles. In Ref. [37] the force distribu¬ 
tion in packings of frictionless particles at jamming was 
found to be singular: 

P(/) oc f 9 t e - f/fo (4) 

as confirmed by several numerical studie^EHMIlll. Note 
that the exponential tail was observed earlier in Ref. [35] 
but is not crucial; a Gaussian cut-off would not change 
our scaling predictions below. The singularity at small 
/ was shown to be necessary for the stabil ity of pack¬ 
ings toward rewiring of their contact networl i 2G F° l More 
precisely, contacts need to be classified in two groups. 
Local contacts have very little mechanical coupling with 
their surroundings, and their density follows Eq.Q. Ex¬ 
tended contacts are coupled to their surroundings, and 
their density follows: 

Pe(/W e =e-//*. (5) 

The latter are less numerous at low force, i.e. 9 e > 9^, 
but more important for our purposd 31 ! Marginal stability 
of packings implie^^l 


2 + 9 e 


( 6 ) 


from which we can compute the effective potential 
Veff({h)) = -1(h)f ((h)) - T\og[(h) + ^(4) 2 + 4T] and 
stiffness 


k((h)) = 


df 

d(h) 


1 (h) 

2 2^/(h) 2 + AT 


(3) 


The effective stiffness is plotted in Fig. UK showing an 
entropic smoothing as T is increased from zero. In the 


in agreement with numerics ind icating 9 1 ss 0.17 and 
9 e ss 0.44. Replica calculations-in d = oo do 
not predict 9t but yield 9 e = 0.42311.P^, consistent with 
the value extracted from numerical observations in finite 
dimensions. 

In what follows, we assume that the distribution of 
Eq.([5]) holds true in the vicinity of the jamming tran¬ 
sition. Note that away from the Maxwell point, P e (f) 
develops a plateau at small forc^HHH. For simplicity, we 
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FIG. 3. Time-averaged gap distribution at T = 10 6 and 
indicated 5(j> (a) as predicted from Eq. §,Eq.Q, and Eq. |8|, 
and (b) as computed from the displacement covariance matrix 
in numerical simulation^ 31 ! 

neglect this plateau in our arguments, since as shown 
previously in the hard sphere regime^ it is not expected 
to affect our results. 

Assuming Eq. 0 wit h 26 ' 291471 9 1 = 0.175 and the ef¬ 
fective potential Eq. 0 immediately yields a prediction 
for the time-averaged gap distribution, as a function of 
/o and T. Using the equation of state derived below, 
this can be written as a function of S(f> 7 plotted in Fig|3|i 
at T = 10~ 6 . As 6<j> is varied, the distribution has a 
complex evolution. Comparing this to numerical results 
from Refj3Tl plotted in Figjsjo, again a very good quali¬ 
tative agreement is obtained, with no fitting parameters. 
Moreover, the scaling behavior of P(h) can easily be ex¬ 
tracted. One finds in particular that P (h) has a max¬ 
imum at W ~ (h), where “ denotes the average over 
contacts, i.e., with respect to Eq.Q. By averaging 
(h) = — / + T/f over contacts, the mean gap is 

W) = -(l + 6 e )fo + T/(f 0 9 e ). (7) 

The height of the maximum of P(/i) scales as P(fi^) ~ 

( fo+T/fo r 1 . 

From P(h) one can also extract the instantaneous co¬ 
ordination Zi ns t by Zi ns t/z = P (h). The result 

is that in most of the phase diagram, Zi ns t < z c : the 
system is instantaneously hypostatic. However, the in¬ 
stantaneous coordination is not the correct variable to 
characterize stability and vibrational properties- this role 
is played by the coordination coarse-grained in time. 

C. Equation of State 

The effective potential and force distribution also im¬ 
ply the equation of state. The Irving-Kirkwood ex¬ 
pression for the stress tensor 55 1 gives its contact contri¬ 
bution as pV/Nc = (1 + h)f, whose thermal average 
is p(V)/Nc = (1 + {h))f. Meanwhile, mechanical sta¬ 
bility of the metastable state allows to write the work 
done in compression from V c in terms of forces and 
gaps. The latter expression, the virtual work princi¬ 
ple derived in Appendix A, equates p(V — V c ) to its 
microscopic expression Nchf ; the thermal average is 



°o 

■^0 °io 0 ° ° o ° ° °o 0 °» 

; :> ° > ° 

1 ° ° ° O*. ° o ° o 

> r<v°c o o\ B v°v^ ; 

°° ° ° o c'° o°o 

\ °"° O 0 000°D° / ' 

\ ^ ' 

5 ° 0 °° o 0 0 ° c 

°^° ° ° o ° ° O °o°o° 

? O \ oo°° o o 0 o 

\A A vV 

: A 


FIG. 4. Illustrative diagram of cutting argument, show¬ 
ing cut bonds in blue (thick lines), (a) Bonds are cut around 
blocks of size LxL, a useful procedure at small T ; (b) When T 
is large enough, the variational argument is improved by cut¬ 
ting instead the fraction q of weakest extended bonds. From 
Ref. El 


p((V) — V c )/Nc = ( h)f . Subtracting these expressions 
we find pV c /Nc = / = (l+0g)f o , using Eq.Q. Similarly, 
we use the effective potential (h) = — f+T/f and Eq. 0 
to evaluate (h)f = —f 2 + T = — (2 + 6 e )( 1 + 9g)fg + T. 
Substituting this into the virtual work principle and sim¬ 
plifying with (V) = V and ( V c — V)/V = Scjj/cj), valid in 
the thermodynamic limit, we find eventually the equation 
of state 

Vc = , = ^ + V(^) 2 + gT(f> 4 
P (1 +9 e )N c /0 2(/> 2 </>(T 1 (2 + 0g) ’ U 

where Ci = 4(2 + ^)/(</> 2 (l + 9g)). Note that / = 
(l + 9t)fo and in all cases p ~ fq. From Eq.Q, we see a 
temperature scale T ~ |<J<?!>| 2 delineating three regimes: a 
hard-sphere like regime (HS), a soft regime (S), and an 
anharmonic regime (AH) as indicated in Table 1. These 
are denoted ‘unjammed scaling,’ ‘jammed scaling,’ and 
‘finite temperature critical’ in Refl39i respectively. 

In fact, we now show that in the soft regime, tempera¬ 
ture plays a role at a smaller , nontrivial scale. The former 
is thus actually two regimes: soft-zero-temperature (SO), 
and soft-entropic (SE). Their scalings are summarized in 
Table H 

II. VARIATIONAL ARGUMENTS 
A. Density of states 

The effective potential defines a harmonic shadow sys¬ 
tem whose density of vibrational states , D(ui), can be 
bound with a variational argument! 14 1 17 l We follow Ref. 
l34l sketching only the salient modifications from the 
T —► 0 arguments presented there. 

The shadow system is equivalent to an elastic network 
with a fixed time-averaged coordination z > z c . Our 
goal is to construct normalized trial displacement fields 
with a small elastic energy E as measured by the stiffness 
matrix of the elastic network. It is customary to use the 
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FIG. 5. Sketch of bound on D(u>) from variational argument, 
(a) for T < T*, and (b) for T > T *, in log-log axes. 


notation E = w 2 , although our predictions, based on 
thermodynamics, apply both to over-damped Brownian 
particles, or to inertial ones. 

If Q orthonormal modes per unit volume can be found 
with a characteristic frequency to, a variational inequal¬ 
ity implies thatP^ D( w) > Q/u>, where here and in the 
remainder of this section we ignore unimportant con¬ 
stants. To construct trial modes, a useful ansatz near 
the Maxwell point is to cut a fraction q <C 1 of bonds, 
creating a density q — Sz/z c of floppy (i.e. zero-energy) 
modes in the cut system. These floppy modes have mo¬ 
tion transverse to all bonds, but motion parallel to bonds 
only at the cut bonds. In the original, uncut system, 
these modes compress or extend springs at the cut bonds, 
but a judicious choice of cut bonds can lead to small en¬ 
ergy and an optimal bound on D(ui). Two ansatze have 
proven useful: one that creates ‘anomalous’ rnodes^, and 
another that creates ‘sparse’ modetP9. First we suppose 
that the springs are at their rest length. 


weakest bonds; we call these modes ‘sparse’. The trial 
modes created by cutting these contacts have a character¬ 
istic stiffness k c , and characteristic force f c , satisfying q = 
J 0 h P(f)df. This leads to f c ~ pq l /( 1+9c \ From Eq.Q 
and Eq.([|) we find k c ~ p 2 q 2 ^ 1+e ‘"> / (T + p 2 q 2 ^ 1+e ^), 
so that the elastic energy of the trial modes is E ~ qk c . 
These two scales are the only ingredients needed from the 
effective potential of the preceding section; in particular, 
the equation of state does not directly affect vibrational 
properties. 

Although in principle both localized and extended con¬ 
tacts could be cut, it was checked in Ref. [M| that 
only the extended contacts give a useful bound for hard 
spheres. We therefore use only them here, taking P(/) 
from Eq.([5]). 

To obtain a useful bound from this ansatz, the char¬ 
acteristic stiffness must be small; in particular, we must 
have T/p 2 3> q 2 /^ l+6 P> (otherwise k c ~ 1). This can 
be written q <C qr with qT ~ (T 1 / 2 /p) 1+9e . In this 
case the modes have a characteristic frequency u(q) ~ 
y/E ~ \Jqk c ~ pg 1 /( 4b )T _1//2 , and the variational in¬ 
equality gives 

D(co) > (T/p 2 ) 2b tu- a = D a {u>), (10) 

with 

a=(l-0 e )/(3 + 0e), 6 = (1 + 0e)/(6 + 29 e), (11) 

as obtained previously for hard sphered^. This holds for 
any q > 8z, implying that ui > ui s with 

~ (p/Vf) Sz 1/4b . (12) 


B. Anomalous modes 


D. Optimal bound 


As discussed in Ref. HI trial modes can be created by 
cutting bonds along blocks of size L ~ 1 /g> 1 (see Fig. 
&)• The induced floppy modes have large motions at the 
cut bonds; by modulating these modes with plane waves 
with nodes at the cut bonds, the resulting trial modes 
have a small frequency oj ~ q/Vk. It turns out that in 
the regimes where anomalous modes can dominate the 
low-frequency spectrum, k ~ 1 and oj ~ q. This leads to 
a bound: 

D(u) > w° = (9) 


At finite T, both ‘anomalous’ and ‘sparse’ bounds are 
useful, in different parts of the phase diagram, and dif¬ 
ferent frequency ranges. The optimal bound on D(ui) is 
simply whichever is larger of Eqs.(|TT){ [?J[) : 

D(ui) = max(£> s (w), D a (u)) (13) 

Accordingly, the characteristic frequency of the lowest- 
frequency modes is oj* = min[w a ,w s ]. The equality oj a = 
uj s defines a characteristic temperature 

T* ~ p 2 Sz a/2b , (14) 


valid above the characteristic frequency uj a ~ Sz. This is 
sketched in Fig(5^,. Since this argument is purely geomet¬ 
rical, it applies without change in the case of an effective 
potential, as considered here. 


C. Sparse modes 


and therefore two cases: 

Case T < T*: the vibrational spectrum is dominated 
by anomalous modes and the density of states is flat, as 
occurs at zero temperature for soft particles. 

Case T > T*: the lowest soft modes are sparse-like. 
However since Eq.( |10[ ) is a decreasing function of w, even¬ 
tually Eq.(10) yields D(oj) > 1, and the bounds coincide. 
This occurs at a frequencjPS 


WT = ( T/p 2 ) 2b / a . 


The second ansatz, discussed in Ref. |3U creates trial 
modes that only extend and compress a fraction q of the 


(15) 
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FIG. 6. Effective medium theory predictions for a soft-sphere 
system at marginal stability, (a) Density of states D(to) at 
T = ICG 8 and indicated values of Srf). The triangle has slope 
2. (b) Shear modulus fi at indicated values of T. 


So the vibrational spectrum contains both features: for 
to s <C to -C wj, D(to) is hard-sphere like and decays as 
a power-law of frequency. However it becomes flat for 
lot < w < 1. We include in the latter case the possibility 
that lot ~ 1, in which case the flat regime is not visible. 
These two behaviors of D(to) are sketched in Fig{5] 


E. Coordination and marginal stability 


The coordination z is a natural variable to describe 
vibrational properties, but it is generally not experimen¬ 
tally accessible. Here we relate it to T and p by consider¬ 
ing that the system is marginally stable, as previously ar¬ 
gued from dynamical consideration^^, and as confirmed 
by replica calculations in infinite dimension^. The elas¬ 
tic energy of a spring network under compression has two 
components: a positive (stabilizing) one, due to motion 
parallel to springs, and a negative (destabilizing) one, 
due to motion transverse to springs. Marginal stability 
implies that for the softest sparse or anomalous modes, 
the destabilizing transverse part of the elastic energy is 
of the order of the stabilizing part. The most unstable 
(anomalous or sparse) trial modes used to bound D(to) 
are those at frequency to*, with energy to* 2 . For such 
modes, the destabilizing energy scales as the pressure— 
—p. We thus get to* ~ ^Jp. This sets Sz and T* as 


5z 


( T/p) 2b t>T* 
y/p T <T* 


(16) 


and 


T* 


2 — a 

P 1 -“- 



FIG. 7. Effective medium theory prediction for 5z(T,5rf>) at 
marginal stability, at various temperatures from T = 10 -1 ° 
(blue) to T = 10 -5 (purple). Here we take 9 e = 0.42. 


F. Mean-squared displacement 


These results yield a bound on the particles’ mean- 
squared displacement ( SR 2 ). 


(SR 2 ) 

T 


J w Jw>u* 


D M 


dto 


> 


(T/p 2 ) 2 V 4b - 2 

1 /to* 


T > T *, 

' j 1 ^ 'j 


> 


(T/p) 1+2b 

T/p 1 / 2 


T > T*, 
T <T* 


(18) 


From effective medium theory, which makes detailed pre¬ 
dictions for the shape of D(to) for to < to* (see below), 
we predict these bounds to be saturated. 


III. PHASE DIAGRAM (0,T) 
A. Soft regime (SO), (SE) 


At T = 0 and (j) > <fi c , the system is in a ‘jammed 
scaling’ regime 39 with p ~ d<j>, Sz ~ ^/p, a nd a plateau 
of ‘smooth’ modes above to* ~ This soft- 

zero-temperature (SO) regime persists for T < T* ~ 
5(j) < S 2 ~ a SR 1 ~ a \ at which point the bound (101 first be¬ 
comes applicable. For T* < T < S(f> 2 , the density of 
states has the form (flO]) , and Sz ~ ( T/p ) 2b , but the pres¬ 


sure is still dominated by contacts. This is the nontrivial 
regime in which typical contacts are elastic, but enough 
weak contacts are entropic to affect vibrational proper¬ 
ties. We call it SE - soft entropic. 

Collecting results, and using Eq.(fl8|), we get 


and 


with 


w* ~ 64 i 1/2 (SO), (SE) 


(SR 2 ) > 


T Sc/- 1 / 2 

(T/s<t>y 


(SO) 

(SE)' 


3 cl 4 + 2 0 e 
2~2 = 3 + <9 e ’ 


(19) 


( 20 ) 


(17) 


K = 


( 21 ) 
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B. Anharmonic regime (AH) 

Both above and below jamming, when the temperature 
is large enough, T > Scjj 2 , the pressure is entropic: p ~ 
\/T. In this regime 

uj* ~ T 1/4 (AH) (22) 



and 

(5R 2 ) > T K/2 (AH) (23) 


FIG. 8. Phase diagram of Fig. [I] in (p, T) space. Lines 
and colors are as in Fig. [l] In these parameters, the regime 
AH folds into the curve p ~ %/T. Note that the lines are 
schematic. 


C. Hard-sphere regime (HS) 


When 8cj> < 0, as T — > 0 the system behaves like hard 
spheres, withp ~ T/\6<f>\, Sz ~ \S(jj\ 2b , and the bound (10) 
above uj*. These scalings persist until the anharmonic 
regime T ~ |50| 2 . Therefore we find simply 


uj* ~ VT\ S(j)\- 1/2 (HS) 

and 

(SR 2 ) > \5(t>\ K (HS) 


(24) 


(25) 


The results (20) and (25) display a striking asymmetry: 
the nontrivial temperature scale T*, although it can be 
defined both above and below jamming, plays no role in 
the latter case. Physically, this is because when 5(j> > 
0, there must be a temperature scale where the softest 
modes transition from ‘anomalous’ to ‘sparse,’ whereas 
when 8 <j) < 0, the softest modes are always ‘sparse’. 


IV. EFFECTIVE MEDIUM 

The scaling predictions of the previous section can 
be derived and extended with effective medium theory 
(EMT), a m ean-fi eld approximation to self-consistently 
treat disordeJ 15 5 3ES. Although EMT is formulated at 
T = 0, an effective potential facilitates its use at any T. 
An arbitrary distribution of disorder is allowed: here we 
assume that forces are distributed as Eq. © , as in the 
variational bounds above. We follow the same proce¬ 
dure as our previous work^IES, discussing all details in 
Appendix B. 

In addition to reproducing the previous results, with 
values for prefactors, EMT also gives the behavior of the 
shear modulus and density of states at all uCl. Here we 
focus on the results for d = 3 at marginal stability, and 
only work to leading order in Sz. By numerical solution 
of the EMT equations at marginal stability, we obtain 
Sz(T,<j>), plotted in Fig. [7j The asymmetry between 
S(f> > 0 and 8(f> < 0 is clearly visible. A curious feature is 
that at fixed T, Sz has a minimum at a finite 5(f>. 

The EMT prediction for D(uj) is plotted in Fig. [6ji. 
Above u*, it behaves as described in the foregoing text. 


Below uj*, we find D(uj) ~ uj 2 ", we note that our previ¬ 
ous worlPS erroneously predicted D(uj) ~ uj 2+a in this 
regime: an incorrect approximation was used to obtain 
scaling behavior below uj* , affecting only this prediction. 

EMT predicts that the zero-frequency shear modulus 
p oc —Sz + \JSz 2 + Cp, for a constant C. It can be 
checked that in all cases p < Sz 2 , so that 

p~p/Sz, (26) 


cons istent with arguments presented earlier in the T —> 0 
limiiPMSHSU, anc j w ith the replica theorjf^. This implies 

Scjj 5 (SO) 

(SE) 

T\S(j)\- K (HS) 1 ’ 

T 1 -^ (AH), 

plotted in Figj6]). In all cases we find 

p(SR 2 ) ~ T. (28) 



We note that this shear modulus only describes relax¬ 
ation within a metastable state. 

Finally, EMT predicts the spatial correlation of mode 
displacements, £ s (uj). 

1 / \JujD(uj), leading to 


L(uj*) 


For uj > 

uj*, we i 

S<j>~i 

(SO) 

( T/Scjj)^ 

1 (SE) 


(HS) 

tV 

(AH). 


(29) 


In principle, t s can be measured by static response, as 
shown previously for T = 0 soft solidiPS. 


V. COMPARISON WITH SIMULATIONS Si 
EXPERIMENTS 

Packings of soft harmonic spheres at finite tem¬ 
perature hav e been studied numerically by several 
groupPi .3 3J36HH1 A major goal of these works, partic¬ 
ularly Refs. I5U1I5T1 andSU was to show that vibrational 

















properties become strongly anharmonic as packings are 
heated. As evidence, it was observed in Refs. [30j[3lj and 
ESSthat as a jammed solid at 6<j> > 0 is heated, the plateau 
in D(uj) develops a negative slope at low frequency. We 
have shown that this phenomenon can be understood 
within a harmonic effective theory, and simply reflects 
the increasing importance of weak contacts (and there¬ 
fore ‘sparse’ modes) as the effective potential develops a 
small-force thermal tail. 

In Refs. I3B, !35, and mi the three simple regimes HS, 
S, AH were observed, with the scaling T ~ 5(j) 2 predicted 
hereP. However, neither of these works found the sub¬ 
tle distinction between SO and SE regimes. Because the 
temperature scale T* ~ \S(/)\ 2 - 20 (using 9 e = 0.42) is close 
to Sc/) 2 , precise measurements of D(uj), (SR 2 ), or /x are 
needed to test our predictions. However our non-trivial 
prediction on the shape of the density of states (which 
can present a plateau or a power-law decay depending on 
the control parameters and the frequency) appears to be 
confirmed numerically, see e.g. Fig.2 in Refl551 and Fig.3 
in Ref. 1551 

Several colloidal experiments have aimed at mea¬ 
suring fi nite-te mperature vibrational properties near 
jammin gj 44 l 45 ^ , although these observations may not 
be close enoug h to the critical point to obtain accurate 
exponents 39 6 “ However, recent experiments on granular 
medicP^ 70; may probe this regime, albeit in the presence 
of friction. In Refs. |55] and [7DJ heterogeneity associated 
with dynamics of contacts was used to infer Widom lines 
emerging from a common point at T = 0, <f> = <j> c , con¬ 
sistent with the results of this work. A more detailed 
comparison could be possible by computing dynamical 
susceptibility, along the lines of Ref. [35J 


VI. CONCLUSION 

We have proposed a description of the vibrational 
properties of thermal soft spheres near the jamming tran¬ 
sition, based on a real-space description of soft vibra¬ 
tional modes. These modes are of two types - sparse and 
anomalous- depending on packing fraction, temperature 
and on the frequency considered. Our work shows that 
even when anharmonicities are extremely strong, they 
can be tamed by time-averaging, allowing one to define 
an effective potential from which static linear response 
and thermal fluctuations can be computed. Ultimately 
all our scaling predictions can be expressed in terms of 
one exponent 6 e characterizing forces in jammed pack¬ 
ings, which can be extracted numerically. 

9 e is also very well estimated by replica calcula¬ 
tions in infinite dimensionIn the hard sphere case 
(cf) < (f> c , T —» 0) such calculations lead to identical re¬ 
sults for th e me an-square displacement^ and the elas¬ 
tic moduluj 65 * 71 ( However this approach currently does 
not fix the exponents for soft sphered^, nor does it ac¬ 
cess the shape of the vibrational spectrum nor the length 
scale characterizing normal modes. It would be very in¬ 


teresting to extend the replica method to these questions, 
enabling a full comparison with the present results. 

Although we have focused on the proximity of the jam¬ 
ming transition, we expect our description based on an 
effective potential to be qualitatively appropriate away 
from it. For example, if the packing fraction is decreased 
(or the temperature incr eased), an elastic instability is 
expected to occui^lISlMI near the glass transition where 
configurations become unstable. Denoting cf >o such a 
packing fraction, effective medium 15 predicts that the 
elastic modulus jumps and displays a square-root singu¬ 
larity at that point fi ~Uo + Cy/(f> — (f> o, a s also found in 
Mode Coupling Theorjf2212II, replica tlieorjUHIll, an d the 
random-first-order-transition scenario' 4 . Furthermore, 
numerics support that an elastic length scale diverges in 
two dimensions a^l l c ^ l/((f>— t^o) 1 ^ 4 and therefore de¬ 
creases as the system enters the glass phase. We expect 
these behaviors to be smoothed to some extent, since in 
the vicinity of the glass transition the separation of time 
scales required to define an effective potential disappears. 
Direct extraction of elastic length scales as a function of 
packing fraction or temperature would be very informa¬ 
tive to test this prediction. 

Finally, one may wonder what happens to the soft 
sphere phase diagram when the system is driven by an 
imposed shear. For T = 0 and (f> < (f> c it has been shown 
that the dynamics of over-damped hard spheres is gov¬ 
erned by an operator whose low-energy modes are of the 
anomalous typ^^, so that velocity correlations decay pri- 
marily on a length scale l c ~ 1/ \Jz c ~ z independent of 
9 e 75 7 ® : . It would be interesting to see if the notion of 
effective potential could be used to extend such consid¬ 
erations to finite temperature, a path recently proposed 
in Ref. [771 , 
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In these Appendices we (A) derive the effective poten¬ 
tial, and (B) describe the effective medium theory. 


Appendix A: Effective Potential 

Here we derive, in more detail, the effective potential 
used in the main text. We consider a system of N soft 
spheres at constant pressure p that is fluctuating around 
an isostatic metastable state. The latter has a Gibbs 
partition function 


Z(p,P) = J DSr 


e -0pV e -pu 


(A.l) 


where U is elastic energy, V is volume, p is pressure, 
f DSf = [I, / dSf, and /3 = 1 /{ksT). The displace¬ 
ments are with respect to a configuration in which all 
particles are just touching, at volume V c . 

Metastability of the state is interpreted to mean that 
although the state may have a finite lifetime, the time- 
averaged positions of the particles are in force balance. 
Then there exist contact forces {f a } balancing the exter¬ 
nal forces {Fi}: 


Fi = ^2f a (n a ), 


(A.2) 


for each i, where (n a ) = (fj — fi)/\(fj — f) | is the time- 
averaged contact vector. Summing over all particles and 
contracting with an arbitrary displacement held, we get 
an expression for the work done in the displacement held: 

W = ^2 Fi ■ Sf = - ^ fa(n a ) ■ (Sfj - Sfi) 

i ot 

= -£/a(fca + 0(/£)), (A.3) 

a 

the virtual work principlfP^. Here we have dehned the 
gap h a between particles, h a = n a ■ (Sfj — Sfi) (with 
h a < 0 for overlap). The error in Eq.(A.3) comes from 
the deviation of time-averaged contact vectors from in¬ 
stantaneous ones; this is small if particle displacements 
are also small, \Sfi\ 1. Because of isostaticity, each 
force f a in this expression can be written f a = pf a , where 
f a is hxed on the mesoscopic time r and depends on the 
time-averaged positions of the particles. The pressure 
simply sets the force scale. 

Isostaticity implies that the number of degrees of free¬ 
dom in the particle displacements {Sf}, dN, is precisely 
equ al to the number of contacts, Nc = zN/2. Thus 
Eq.(A.l), originally an integral over the displacements 
{Sfi}, can be written as an integral over the gaps {h a } 
between contacting particles, i.e., 


Z(p,P) = J DhJ 


e -0pV e -0U 


(A.4) 


where the Jacobian J = | const ant +G(h)\: to leading or¬ 
der, J depends only on the time-averaged positions of the 


particles, with a small correction due to contact vector 
rotation. Since compressing or dilating the system from 
volume V c to V req uires work W = —p(V — V c ), we can 
substitute Eq.(A.3l into Z to obtain 


Z(p,/3) = Je~ 0pVc JJ f dh a e~ p Fho, e - 0 U(h a )^ 


where we assume that U is a pair potential, thus depend¬ 
ing only on h a ■ Strictly, each gap h a should be bounded 
from above by some nontrivial function of the positions of 
the rest of the particles. However, since we are only inter¬ 
ested in small displacements, we ignore this fact and let 
the upper limit on each h a —> oo: large gaps are strongly 
suppressed by the Boltzmann factor. ( A posteriori, one 
can check the smallness of temperature needed for gaps 
to remain C 1.) Finally, we are led to a single-gap par¬ 
tition functioi ^ 32 * 33 ! 


Z(/3,f)=e-P G = j dh e -P u W e -P fh , 


(A. 6 ) 


where the force / fixes the time-averaged gap (h) through 
(h) = dG/df. The effective potential 14// should be a 
function of the time-averaged gap (h); we define it by 
Legendre transform V e ff((h)) = G(f) — f(h), from which 
it follows that / = —dV e ff({h))/d(h), as desired. We con¬ 
sider a finite-range harmonic potential U{h ) = ^e\h/a \ 2 
when —\a<h< 0, and 0 otherwise. We now take units 
in which e = a = ks = 1 . 


Eq.(A. 6 ) and / = —dV e ff{{h))/d{h) define implicitly 
the effective potential. For our choice of U(h), we obtain 
an error function: 


z = f dhe -^ + J_ 


Pf VP 


0 

- 1/2 

!>/■ j je ‘ 


dhe-Pfbe-iPW 2 (A.7) 


dz e 2 " 


'vW-f) 


(A. 8 ) 


We are interested in /3 —> oo so yffiif — 5 ) ~ — 00 . Then 
the equation defining the effective potential is 


m = - 


dlnZ 


df 


(A.9) 


-1 

IT 


1 1 2 r y S~@f 1 

'472 + 1 + /3 1 / 2 /e2 /S/ ' / dze~ 2 
PJ J -00 


We see that the behavior of Eq.(A.9) depends on the 
scaling variable VPf ■ When /3 / 2 1, the contact term 

dominates, while when /3 / 2 <C 1 , the entropic term dom¬ 
inates. The integral is always 0(1), so we approximate 
it by Vn/ 2 . 

Entropic regime: Let /3f 2 <C 1. Then Z ss l/(/3/) 
and 


(h) 


1 

Jf 


(A.10) 
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This recovers the hard-sphere effective potential derived 
in Refl33l It holds when / <C 1/-//4 or (h) 1/-//L 

Contact regime: Let /3f 2 1. Then Z ss 

\/tt/P/2 and 

(h) «-/ (A.ll) 

This recovers the original harmonic potential. It holds 
when / » l/\//J, or ( h ) <C — l/y/]3. 

Therefore there is a characteristic scale 1 /y/]3 = \/T 
both for gaps and overlaps. We cannot analytically solve 
the equations in the nontrivial regime — \JT < (h) < vT 
but there is a smooth cross-over (for example, there is no 
/ with d(h)/df = 0, since d(h)/df oc — ((h — (h )) 2 ) < 0). 

In the main text, we want to use 14// to obtain scal¬ 
ing behavior of vibrational properties; we can do so by 
replacing the cumbersome implicit expression with a sim¬ 
ple form that incorporates the correct limiting behavior. 
The simplest is the one satisfying (h) = —f + T/f, which 
leads to an effective force law 


where P = z/zq. The bond stiffness distribution is 
P(fc) = df /dk P(/). In EMT, the elastic behavior of a 
random material is modelled by a regular lattice with ef¬ 
fective frequency-dependent stiffnesses; here we will have 
a longitudinal stiffness, fell, and a transverse stiffness 
— k ± . Writing ~ for disorder average, the EMT equa¬ 
tions are 


0 1 - (fell - fc a )Gll 1 + (fc- L -/a)G- L ’ (B ' 4) 

where GH and G 2 - are related to the Green’s function 
G(w) = (Ad — mu 2 ) 1 by 


= fi a ■ (a|G|o!) • n a 


G = 


d- 1 L 


tr((a|G|a)) - G 11 


(B.5) 

(B.6) 


with (a| = (z| - (j\. 
EMT equations 


Using Eq.(A.12) one obtains the 


/«/>» = -^) + 5V / W 2 +4T, (A.12) 

from which we can compute the effective potential 
Veff((h)) = -1(h)f ((h)) - Tlog[(h) + vWT4T] and 


stiffness 


k((h)) ee 


df 

d(h) 


1 (h) 

2 2 a/ (ft ) 2 + 4T 


(A.13) 


o= +p 


1 — fell Gil 


1 + 


CfTfl 


-1 | 1 

Gil + Gil 2 (1 + c) 

-[df ^ 


cTf 0 " 2 + P 


(l-P)fc 1 - 

0 = ----h P 

1 + fcJ-G-L 


1 

G 1 


Cffo 
G x2 


df 


fe 


-f 


C2 ~ f _ 


These relations are valid for T -C / 2 and T f 2 , but 
perhaps not close enough to T ~ f 2 . 


Appendix B: Effective Medium Theory 


Effective medium theory (EMT) is a self-consistent ap¬ 
proximation scheme to treat disorder. Here we follow our 
previous work^LitHJ, th e only difference with respect to 
Ref. [34] be ing th e effect of softness on the effective po¬ 
tential Eq.(A.12) (in that work we had f((h)) = T/(h)). 
Accordingly, we only sketch the derivation of our results. 

Equations: We assume that contact forces have a dis¬ 
tribution 


P(/) = C f f 8 e -f'f°, 


(B.l) 


and that forces and bond stiffnesses follow the effective 


potential Eq.(A.12), Eq.(A.13|. We fix the bond length 


to one, since incorporating the dependence of bond length 
with force leads to subdominant corrections near jam¬ 
ming. A random elastic network of coordination 2 is 
constructed by random dilution of a regular lattice of co¬ 
ordination zo down to z. The stiffness in contact a, k a , 
and the force in the contact, f a are random variables 
distributed according to 


(B.2) 

(B.3) 


where c = (1 — feHGH)/GH, c = c/(l + c), and C 2 = (1 + 
k- L G J -)/(f 0 G J -). Near the Maxwell point we expect |c| <C 
1 and C 2 1 (which can be checked a posteriori ), leading 
to 


o = t iic»-p + < 1 - t '', G "> f 


Gil 


1 + 


„ T ,_ 1 _ e /T(l-fcllGll)Y 
ClTfo —Gii— ) 


0 = k ± G ± - P 


foG ± (d + 1) 
l + k^G 2 - 


(B.7) 

(B.8) 


with a = (6 — l)/2, Gi = 7 r/( 2 T 1+ g sin | 7 ra|), T t = 
/ 0 °° x t ~ 1 e~ x dx. 

For the Green’s function we consider 



d d q e^ 

p( 27r) d (fell — kpq 2 — mui 2 ’ 


(B.9) 


where BZ = {q : |</| < A} is an approximate first Bril- 
louin zone, 6 is the identity tensor, p = N/V is the num¬ 
ber density, and k 2 - = (d — 1 )k 2 ~. Isotropy of G implies 
an identity 


G 11 = G 2 - 


2d 1 
zo fell — k ± 



mur 

d 


tr(G(0, ui)) . 


P EMr(k a ) = (1 - P)S(ka) + P P (k a ) 
P EMT(P) =(l- P)6(f a ) + P P(/a) 


(B.10) 
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Assuming w <C yJkW jm and d > 3, it can be checked that 




(B.ll) 


with 


Ai = 


2n d/2 A d-2 

pr d/2 (27r) d d- 2 


(B. 12) 


Eqns.(B.7),(B.8l,(B.10),(B.ll) are so i ve d f or § z <g; i by 
supposing 


fell ~ Sz^, k ± ~ Sz v , (B.13) 

ui ~ Sz**, T ~ Sz", (B.14) 


this case, leading to FigjTJin the main text. Since fco = 
fell (oJo ) 5 at marginal stability fco — fcA = fco(l + 0(Sz)) is 
precisely the shear modulus /i. Hence 


A* = 


7i 

272 


Sz+ yj Sz 2 + 73/0 , 


(B.20) 


as discussed in the main text, and plotted in FigjGjo. 
Density of states: The density of states 


D(uj) = ( 2 w/ 7 r)Im[tr[G( 0 , w)]] 

= 2 cL 4 ^Im[l/Afc] + ... 
7 r 

2dA\u: Im[Afc] 

7T | A A’ | 2 + 


(B.21) 


and balancing terms in the above equations. This nec¬ 
essarily leads to a description of the behavior near T*, 
where both hard-sphere and soft-sphere physics is present 
at small frequency. We find 


V = 7 = = 2C = 2, v = , (B.15) 


recovering the scalings described in the main text. 
Results: The leading order equations for AA,fcH are 


*A = Pf 0 (l+9)~p (B.16) 


and 


0 = z c Aiuj 2 — k^Sz + k^z c + ak + c 3 fc^ 


/TfcH \ 


a+1 


V /o 2 ) 


(B.17) 


with a = z 0 - z c and c 3 = Ci 2 0 (l - z c /z 0 ) a+1 (z 0 /z c ) a . 
This equation contains the leading-order terms (in 8z) 
for both the harmonic T = 0 theory from Ref. [T5] and 
the hard-sphere theory from Ref. [34] When all terms are 
important, it describes behavior near T*. Other regimes 
are obtained when one or several terms becomes small. 
We cannot analytically solve Eq.(|B.17|), but we can de¬ 


termine some of its key properties. 

For example, we expect A (uj) to have an onset fre¬ 
quency u >0 where the density of states rises from zero. At 
this point we must have \dk^/du\ = 00 , implying 


0 = 5z — 2 oAq — (a + 2 )c 3 


Tk[ J 
. fo 


a +1 


(B. 18 ) 


and 


0 = 72*0 


+ 71 5zk\ - z c (AiWp + /A), 


(B.19) 


where k 0 = fcH(u; 0 ), with 73 = (1 + 9 )/(3 + 6), 72 = 
a(l — 8 )/(3 + 9 ). fcp a nd ujq can be found by numerical 
solution of Eqs.( |B.18| ),(B.19) for a prescribed value of Sz. 

Marginal stability: The condition of marginal stability 
corresponds tcP-^ ujq = 0. Therefore 8 z(p,T) at marginal 
stability can be found by solving Eqs.(B.18),(B.19) in 


with Ak = A — fcA. We determine it by numerically 
solving Eq.(B.17) at each oj. It is plotted in Figjh^. 

Scattering Length: Finally, we can extract the asymp¬ 
totic behavior of the Green’s function for large r. To 
leading order, log(G(r, w)) ~ —r / t s (oj) +iur / where 

= — uW 1 1Afc|/Im[\/Afc] and v(u>) = |Afc|/Re[VAfc] 
are, respectively, the scattering length and sound velocity 
at frequency u. 

To extract t 8 at the characteristic frequency w*, we 
use our knowledge o f th e shape of Afc ~ k" known from 
our previous works 1 ®^, which is confirmed by numeri¬ 
cal solutions to Eq.(B.17): near 0 A, we have Re[feH] ~ 
constant ~ p, while — ImffeU] >Re[feH]. This implies that 
|Afc| ~-Im[feH]. Using 2Im 2 \/Afc = |Afc|—Re[Afc] ss |Afc| 
along with the definitions of i s , and Eq.(|B.21 ), we find 


4(<A) = - 


I Afc | 


w*Im[\/Afc] 

lAfcl 1 / 2 


P(w> 1 


- 1/2 


(B.22) 


from which we can determine the scaling properties of 
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